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ABSTRACT 

The detection of the primordial _B-mode spectrum of the polarized cosmic microwave background 
(CMB) signal may provide a probe of inflation. However, observation of such a faint signal requires 
excellent control of systematic errors. Interferometry proves to be a promising approach for over- 
coming such a challenge. In this paper we present a complete simulation pipeline of interferometric 
observations of CMB polarization, including systematic errors. We employ two different methods 
for obtaining the power spectra from mock data produced by simulated observations: the maximum 
likelihood method and the method of Gibbs sampling. We show that the results from both methods 
are consistent with each other, as well as, within a factor of 6, with analytical estimates. Several 
categories of systematic errors are considered: instrumental errors, consisting of antenna gain and an- 
tenna coupling errors, and beam errors, consisting of antenna pointing errors, beam cross-polarization 
and beam shape (and size) errors. In order to recover the tensor-to-scalar ratio, r, within a 10% 
tolerance level, which ensures the experiment is sensitive enough to detect the B-signal at r = 0.01 in 
the multipole range 28 < ^ < 384, we find that, for a QUBIC-like experiment, Gaussian-distributed 
systematic errors must be controlled with precisions of |(7rms| = 0.1 for antenna gain, |erms| = 5 x 10~^ 
for antenna coupling, 5rms ~ 0.7° for pointing, C,rrns ~ 0.7° for beam shape, and ^.rms = 5 x 10~* for 
beam cross-polarization. 

Subject headings: cosmic background radiation - cosmology: observations - instrumentation:interferometers - methods: 

data analysis - techniques: polarimetric 



1. INTRODUCTION 

The Cosmic Microwave Background (CMB) has be- 
come one of the most fundamental tools for cosmology. 
High-precision measurements of the CMB polarization, 
especially d etecting the primordial "B -mode" polariza- 
tion signals ([Kamionkowski et al.lll997l) . will represent a 
major step towards understanding the extremely early 
universe. These S-modes are generated by primordial 
gravitational waves. A detection of these signals would 
probe the epoch of inflation and place a n important con- 
strain t on the infiationary energy scale (jHu fc DodelsonI 
|2002|) . In addition, the secondary B- modes induced by 
gravitational lensing encode information about the dis- 
tribution of dark matter. However, the i3-mode signals 
are expected to be extremely smal l and current experi- 
ments can only place upper limits (jHinshaw et al.|[2012D 
on the tensor-to-scalar ratio; the quest for the B-modes 
is a tremendous experimental challenge. 

Due to the weakness of the i3-mode signals - the largest 
signal of the primordial _B-modes is predicted to be less 
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than 0.1/xK - exquisite systematic error control is cru- 
cial for detecting and characterizing them. Compared 
to imaging systems, interferometers offer certain advan- 
tages for controlling systematic effects because: (1) an 
interferometer does not require rapid chopping and scan- 
ning (Timbie et al. 2006) and, with simple optics, inter- 
ferometric beam patterns have extremely low sidelobes 
and can be well understood; (2) interferometers are in- 
sensitive to any uniform sky brightness or fluctuations in 
atmospheric emissions on scales larger than the beam 
width; (3) without differencing the signal from sepa- 
rate detectors, interferometers measure the Stokes pa- 
rameters directly and inherentl y avoid the l eakage from 
temperature into polarization (jBunnl I2007D caused by 
mismatched beams and pointing errors, which are se- 
rious probl ems for i?-mod e detection with imaging ex- 
periments (Hu et al. 2003; Su et al. 2011; Miller et aLl 
[^09; O 'Dea et al. 2 007; Shimon et al. 2008; Yadav et al] 
2010: ,Takahashi et al.ll2010f ): (4) for observations of small 
patches of sky, E-B mode separation would be cleaner in 
the Fourier domain for interferometric data than in real- 
space; and (5) with the use of redundant baselines, sys- 
tematic errors can be averaged out. In addition, they of- 
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fer a straightforward way to determine the angular power 
spectrum since the output of an interferometer is the vis- 
ibility, that is, the Fourier transform of the sky intensity 
weighting by the response of the antennas. 

Interferometers have proved to be powerful tools 
for studying the CMB temperatur e and polar i zation 
power spectra. In fact, DASI (jKovac et al.l |2002|) 
was the first instrument to detect the faint CMB po- 
larization anisotropics. Pioneering attempts to mea- 
sure the CMB temperature anisotropy with intcrfcr- 
ometers were m ade in the 1980s (Martin ct al. 1980; 
Fomalont et al. 1984; Knok e et al.lll984l : iPartridge et aO 
1988t iTimbiefc Wilkinson 1 19881) . Several groups have 
successfully detected the CMB anisotropics. The CAT 
telescope was the first interferome ter to actual l y de- 



tect structures in the CM B (O'S uUivan et al.l Il995l: 



Scott et al.lll996l : lBaker et al . 1999). CBI (Pearso n et al 
2003^ and VSA (jPickinsoii et al. 2004; Gr ainge et al 
20031 ) have detected the CMB temperature and po- 



larization angular power spectra down to sub-degree 
scales. In the next few years, the QUBIC instru- 
ment (jOubic Collaboration et aLl 120111 ) based on the 
novel concept of bolometric interferometry is expected 
to constrain the tensor-to-scalar ratio to 0.01 at the 90% 
confidence level, with 1-year of observing. 

On the theory side, the formalism for analyz- 
ing interfer ometric CMB da ta has be en well- 
devel op ed (iWhite et al.l 119991: iHobson fc M aisinger 



2()0l ^IPark et al.l 120031: IMvers et al.l 1200.1 l26oft 
Hobson fc Magueijol Il996( ). A pioneering study of 
systematic effects for interferometers based on_an 
analytic approach has been performed bv iBunnI (|2007[ ). 
However, this approach is of course only a first-order 
approximation for assessing systematics, since many 
important effects, such as the configuration of the array, 
instrumental noise, and the sampling variance due to 
finite sky coverage and incomplete uu-coverage, are not 
taken into account. Any actual experiment therefore 
will naturally require a complete simulation to assess 
exactly how systematic effects bias the power spe ctrum 
recovery. In this concern, iZhang et al.l ()20 12f) have 
presented a simulation pipeline to assess the systematic 
errors, mainly focusing on pointing errors. With a full 
maximum likelihood (ML) analysis of mock data, the 
simulation agrees with the analytical estimates and 
finds that, for QUBIC-like experiments, the Gaussian- 
distributed pointing errors have to be controlled to the 
sub-degree level to avoid contaminating the primordial 
i?-modes with r < 0.01. 

Nevertheless, a comprehensive and complete analysis 
of various systematic errors on CMB power spectrum 
measurements has not been undertaken so far. In this pa- 
per, therefore, we perform a detailed study to completely 
diagnose the most serious systematic effects including 
gain errors, cross-talk, cross-polarization, beam shape 
errors and pointing errors, on the entire set of CMB 
temperature and polarization power spectra. In order 
to assess the effect of the systematic errors on _B-mode 
detection and set allowable tolerance levels for those er- 
rors, we perform simulations for a specific interferomet- 
ric observation with an antenna configuration similar to 
the QUBIC instrument . We also extend the analytical 
expressions ()Bunnll2007[ ) for characterizing systematic ef- 
fects on the full CMB power spectra. 



For verifying the power spectrum analysis, we employ 
two completely independent codes based on the Gibbs 
sampling algorithm and the maximum likelihood tech- 
nique. The use of Gibbs-sampling based Bayesian in- 
ference to interferometric CMB ob servatio n s has been 
successfully d e mons trated by iSutter et al.l ()2012f ) and 
iKarakci et al.l (|2013f ). It allows extraction of the under- 
lying CMB power spectra and reconstruction of the pure 
CMB signals simultaneously, with a much lower compu- 
tational complexity in contrast to the traditional maxi- 
mum likelihood technique (Hobson & Maisinger 2002). 

In this paper, for given input CMB angular power 
spectra, we simulate the observed Stokes visibilities in 
the flat-sky approximation. We believe that the flat-sky 
simulations are sufficiently accurate for the study of sys- 
tematic errors. First, in our simulation, we assume single 
pointing observations with 5° beam width, correspond- 
ing to a sky coverage fraction of fsky = 0.37%. This 
sky patch is small enough to permit the use of the flat- 
sky approximation. Second, all the data analysis pro- 
cesses are established using the flat-sky approximation 
while the mock visibility data are also simulated using 
this approximation. Therefore, a self-consistent analy- 
sis is performed. However, when using a patch cut from 
the projection of spherical sky onto a flat image as an 
"input" map, one s hould take into account the contami- 
nation ()Bunnll2011l) of "ambi guous" modes arising from 
incomplete sky coverage and thus requires an appropriate 
data analysis method to apply for this situation. 

This paper is organized as follows. In Sec. 2, we briefly 
summarize the effects of a variety of systematic errors on 
interferometric CMB observations and describe the an- 
alytical method for estimating those errors. In Sec. 3, 
we describe the simulations interferometric visibilities 
that include systematic errors. In Sec. 4, we review the 
data analysis methods used in this paper, including the 
Gibbs sampling technique and the maximum likelihood 
approach. In Sec. 5, we assess the systematic effects on 
the CMB power spectra. Finally, a discussion and sum- 
mary are given in Sec. 6. The appendix contains the 
complete analytical expressions for the systematic effects 
on the full CMB power spectra. 

2. SYSTEMATICS 
2.1. Instrument Errors and Beam Errors 

In a polarimetric experiment, the Stokes parameters 
I,Q,U and V can be obtained by using either linear or 
circular polarizers. For a given baseline Ujk = — x-,-, 
Xfc being the position vector of the k*^ antenna, t he vis- 
ibilities can be written as a 2 x 2 matrix V^fe (jBunnI 
t2007(); 



V,k = J (fv Afc(f) R • S • R-iA](f) 6-2"".*.-, (1) 
where the 2x2 matrix A^^ (f ) is the antenna pattern and 
I + Q U + iV 



S = 



U~iV I-Q 



(2) 



For a linear experiment, R is the identity matrix and for 
a circular experiment, 
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(3) 



Various systematic errors can be modeled in the defi- 
nition of tlie antenna pattern as follows (Bunn 2007) 



Afc(f) = Jfc•R•A^(f)•R- 



(4) 



where the Jones matrix 3k represents the instrumental 
errors, such as gain errors and antenna couplings. The 
matrix is the antenna pattern that models the beam 
errors, such as pointing errors, beam shape errors and 
cross-polarization. In an ideal experiment = I, where 
I is the identity matrix, and the antenna pattern is given 
as Ag(f) = A(f)I, where Air) a circular Gaussian func- 
tion. 

In this paper we will consider only two types of in- 
strumental errors; antenna gain, parametrized by and 
and couplings, parametrized by e\ and 62 - The cou- 
pling errors are caused by mixing of the two orthogo- 
nally polarized signals in the system. To account for the 
phase delays, the parameters g and e are given as com- 
plex numbers. The Jones matrix for the fc*'' antenna can 
be written as (Bunn 2007) 
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(5) 



For the beam errors, we will consider that each an- 
tenna has a slightly different beam width, ellipticity 
(beam shape errors) and beam center (pointing errors), 
as well as a cross-polar antenna response describ ed by 
off-di agonal entries in the antenna pattern matrix (jBunnl 
|2003); 



A^=4(p,, 



1 

1 + oMfcf^ cos 20 



sin 20 



jM/cfr sin2(?!) 1 - i/ifcfr cos 20 



(6) 

where {p, 0) is an elliptical Gaussian function writ- 
ten in polar coordinates {p,(l)), o is the width of the 
ideal beam and pk is the cross-polarization parameter 
of the /c"* antenna. This particular form of the cross- 
polarization occurs, with p^ = o'^/2, when the curved 
sky patch is projected onto a plane. 

2.2. Control Levels 

The effect of errors on the power spectra can be de- 
scribed by the root-mean-square difference between the 
actual spectrum, C^J^^j, which is recovered from the 
data of an experiment with systematic errors, and the 
ideal spectrum, C^^J^j , which would have been recovered 
from the data of an experiment with no systematic er- 
rors; 
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XY \ 
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(7) 



where X,Y ^ {T,E,B}. 

The strength of the effect of systematics can be 
qua ntified by a tole r ance pararneter a^^ defined 
bv (lO'Dea et al.l [2007t iMiller et all [200I iZhang et al.l 
[2012) 
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where ct^^ is the statistical l-cr error in XF-spectrum 
of the ideal experiment with no systematic errors. 

The main interest in a _B-mode experiment is 
the tensor-to-scalar ratio r which can be estimated 
as (jO'Dea et al.|[2007l ) 



Eft r<BB(r<BB r<BB \/(„BB \2 



(9) 



where b denotes the power band, Cj^^^.^^ is the _B-mode 

spectrum due to weak gravitational lensing and de- 
pends linearly on r through the amplitude of the primor- 
dial -B-mode s. The tolerance pa rameter of r is given by 
ar = Ar/ur (|0'Dea et al.ll2007D ; 



Ar 



-1/2 



(10a) 



(10b) 



For good control of systematics, the value of ar is re- 
quired to stay below a determined tolerance limit. 

2.3. Analytical Estimations 

Analytical estimations of the effect of systematic errors 
on the pol ariza t ion po wer spectra are extensively exam- 
ined in ,Bunnl ()2007[ ). Defining a vector of visibilities 
V — {Vi,Vq,Vu) corresponding to a single baseline u 
pointing in the x direction, for an ideal experiment, we 
can write 

(Ua) 
(lib) 
(He) 
(lid) 
(lie) 
(llf) 

where c^, and c are averages of cos^(20), sin^(20) and 
cos(20) over the beam patterns: 

/|2S(k-27ru)Psin^(20)d^k _^ - 

/|^2(k)|2d2k 

where A^ is the Fourier transform of the ideal beam 
pattern squared. The unbiased estimator for C'^^ — 

C^^ ) is obtained as 
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(7^^ = vt . N 



XY ■ V 



(13) 



where Nxy is a 3 x 3 matrix involving and c (see Ap- 
pendix) . For a baseline pointing in an arbitrary direction 
the analysis is done in a rotated coordinate system: 
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v^ot = cos 26* sin 26* v, (14) 
V -sin26' COS26I / 

9 being the angle between u and the a:-axis. 

The effect of errors on visibihties can be described, to 
first order, by actual = ^ ideal + 5^- Combining ^^ ideal 
and (5v into a 6-dimensional vector w = (v, Jv), we can 
write the first order approximation as (|Bunnll2007| ) 

{^C^Z,)^ =Tr[{UxY-M^f] + {Tr[UxY-Mu.])\ (15) 
where M.w — (w ■ w^) is the covariance matrix of w and 

The error on a particular band power is, then, given as 
an expansion in terms of ideal power spectra: 

{^^rms,bf — Prrns ^ I^XYJ ,J^b (l^) 
I, J 

where p is the parameter that characterizes the error, 
such as gain, g, coupling, e, or cross-polarization, /i, and 
J, J = {TT, TE, EE, BB}. This expression is valid for a 
single baseline. For a system with nf, baselines in band 6, 
AC,^ b must be normalized by 1/ ^/rlZ, assuming there is 
no correlation between error parameters of different base- 
lines. Analytical estimations of the coefficients Kxy i j 
for various systematic errors are presented in the Ap- 
pendix. 

3. SIMULATIONS 

The input /, Q and U maps are constructed over 
30-degree sq uare patches w ith 6 4 pixels per side as 
described in iKarakci et al.l ()2013D with the cosmolog- 
ical par ameters consistent with the 7-yea,r resu lts of 
WMAP f Larson et ani2011l : iKomatsu et al.ll2011| ). The 
tensor-to-scalar ratio is taken to be r = 0.01. The angu- 
lar resolution of the signal maps is 28 arcminutes, cor- 
responding to a maximum available multipole of imax = 
384. The ideal primary beam pattern, A{r), is modeled 
as a Gaussian with peak value of unity and standard de- 
viation of (T = 5°, which drops to the value of 10^^ at 
the edges of the patch, reducing the edge-effects caused 
by the periodic boundary conditions of the fast Fourier 
transformations. Although the patch size is too large 
for the fiat-sky approximation, the width of the primary 
beam is small enough to employ the approximation. 

The interferometer configuration is a close-packed 
square array of 400 antennas with diameters of 7.89A. 
The observation frequency is 150 GHz with a 10-GHz 
bandwi dth. This configuration is simila r to the QUBIC 
design (jQubic Collaboration et all 120 111) . With this fre- 
quency and antenna radius, the minimum available mul- 
tipole is £min = 28. The baselines are uniformly rotated 
in the M?;-plane over a period of 12 hours while observing 
the same sky patch. The resulting interferometer pattern 
is shown in Figure 1. 

The noise at each pixel for the temperature data is ob- 
tained from the total observation time that all baselines 
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Figure 1. Interferometer pattern created over an observation pe- 
riod of 12 hours by 20 X 20 close-packed array of antennas of radius 
7.89A. 

spend in the pi xel. The noise cova riance for the baseline 
Ukj IS given as (IWhite et al.l[T999l) 

where T^ys is the system temperature, A is the observa- 
tion wavelength, rj^ is the aperture efficiency, is the 
bandwidth, n is the number of baselines with the same 
baseline vector, and ta is the integration time. The noise 
value is normalized by a constant to have an rms noise 
level of 0.015/iif per visibility, yielding an average overall 
signal-to-noise ratio of about 5 for the Q and U maps. 

The systematic errors are introduced by calculating 
the visibilities in each pixel according to Eq. [TJ Each 
antenna has random error parameters for gain, cou- 
pling, pointing, beam shape, and cross-polarization er- 
rors drawn from Gaussian distributions with rms values 

of \grmA = 0.1, \erms\ = 5 X lO'^, Srms = 0.1(7 « 0.7°, 

(rms — O.lcr « 0.7°, and ^rms — 5 X 10^**, respectively. 
Here S is the offset of the beam centers of the anten- 
nas and is the deviation in the beam width along the 
principal axes of the elliptical beams. As the baseline ro- 
tates, the beam patterns of the corresponding antennas 
get rotated as well. Whenever a baseline crosses a new 
pixel, the visibility within the pixel, given by Eq. [1] is 
calculated again with the rotated beam patterns. The 
data in a given pixel is taken as the average of all the 
visibilities calculated in that pixel. 

In a circular experiment, the Stokes variables Q and 
U can be simultaneously obtained for the same baseline. 
Thus, for a circular experiment, p^j^c ~ P^irc- However, 
for a linear experiment, direct measurement of Q requires 
perfect cancellation of the much larger / contribution in 
Eq. El Practically, a linear experiment only measures U. 
Since U ^ Q under a 45° clockwise rotation, Q can be 
measured by measuring U with 45°-rotated linear polar- 
izers. Since Q and U are not measured simultaneously by 
the same baseline, in general, the error parameters pf^^ 
and pYin are treated as the distinct parameters in a linear 
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experiment, i.e., pj^^^ ^ ^in- simulate this, we calcu- 
late Vu with a set of error parameters, Then Q and 
U in Eq. [2]are replaced by — i7 and Qj respectively, and 
Vu is calculated again with a different set of parameters, 

"Piin-i to obtain Vq. The simulation requires 4.5 CPU- 
hours for the circular experiment and 13.5 CPU-hours 
for the linear experiment. 



4. ANALYSIS METHODS 



4.1. 



Maximum Likelihood Analysis 

The scheme for the maximum likelihood (ML) analy- 
sis of CMB power spectra from interfcromctric visibil- 
ity m ea surements is pres e nted in [Ho bson & Maisingej 
(pool : IPark et al.l ([20031) : I Zhang et all (|2012D . which 
we briefly summarize here. The ML esti mator of the 
power spectrum has many desirable features (iBond et al.l 
llOQSl'Stuart fc Or Jl987n an d has be en widely applied in 
CMB cosmology ( Bond et al . 1998; Bunn fc Whiti[T997l: 
[Hobson fc MaisingCTn2002[ r 

In practice, we divide the total i!-range into N^, 
spectral bands, each of bin-width A£. The power 
spectrum Ci thus can be parametrized as fiat band- 



powers Cb(b 



1, 



over A£ to evaluate the like- 
[Bond et al.l[T99l 

. In each of the 
to be a constant 
Ai and 



over 



lihood functi on (Bunn & White 199' 
iGorski et all [1996: White et al. 199! 
band-powers, we assume £{£ + l)Ci 
value to characterize the averaged C, 
has Ch = 27r|uf,P5'(|u ft|) as the flat-sky approxima- 
tion ([White et al.l ll999V 

In our case, the CMB signals and the instrumental 
noise are assumed to be Gaussian random fields. There- 
fore, for a given set of CMB band-power parameters 

the signal covariance 



lEE pBB pTE pTB 



matrices can be written as 
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where we introduced the so-called window functions 



'XY given by 



Wz'z'xyH^D 



2tt 



d(f>w wzxt^Z'y^(uj-w)ii*(uj-w) . 



(20) 

where Z,Z' = {I,Q,U} and X,Y = {T,E,B} with 
UJIT = 1, i^UE = sin2(^w, wys = C0s2(^w, i^QE = 
cos20w, ^QB — — sin2(/)w and otherwise zero. 

Due to the fact that the window functions 
^zz'xyH'^I) independent of Cb, the integrals 

of the window functions over w in Eq. only have 
to be calculated once before evaluating the covariance 
matrices. Additionally, if the primary beam pattern 
A{x) is Gaussian, the wind ow functions can be expresse d 
analytically (see d etails in [H obson fc Maisingeil ([2002D : 
IPark et all ( [200l : [Zhanget al. (20l|)). 

We evaluate the likelihood function by varying the 
CMB b and-powers u sing the a bove parametrizati o n. Fo l- 
lowing iHobson fc M aisingc rl (120021): [Park et al.l ([20030 : 
iMvers et all ([2003[ ): IZhang et al.l ([2012f) . the logarithm of 
the likelihood function for interferometric observations is 



given by 

lnC{{C'b}) = nlog7r~log\Cv+CN\-dl{Cv+CN)-^dv , 

. (21) 

where Cy is the predicted signal covariance matrix and 
Cat is the instrumental noise covariance matrix, dy is 
the observed visibility data vector constructed by dy = 
{■ ■ ■ ■,Vi{ui),VQ{ui),Vu{ui);- ■ ■){i = l,...,n) where i 
denotes the visibility data contributed from the pure 
CMB signals and the instrument noise at the i-th pixel 
in the uv-plane an d we have a total of n data p oints. 

As mentioned by IHobson fc Maisinged ([2002( 1 and ref- 
erences therein, the combination of the sparse matrix 
conjugate- gradient technique and Powell's directional- set 
method give a sophisticated and optimized numerical al- 
gorithm for maximizing the likelihood function to find 
the "best-fitted" CMB power spectrum. With an appro- 
priate initial guess to start iteration, independent line- 
maximization is performed for each band-power param- 
eter in turn, while fixing the others. Typically, this pro- 
cess requires a few iterations, of order N^, to achieve the 
maximum-likelihood solution. For about 4000 visibilities 
in a QUBIC-like observation, the maximum-likelihood 
solution of 6 X 6 CMB band-powers can be obtained in 
around 20 CPU- hours. 

Assuming the likelihood function near its peak can be 
well-approximated by a Gaussian, the confidence level of 
the derived maximum-likelihood CMB power spectrum 
is given by the inverse of the curvature (or Hessian) ma- 
trix at the peak. The Hessian matrix is the matrix of 
second derivatives of the log-likelihood function with re- 
spect to the parameters. This matrix is easily evaluated 
numerically by performing second differences along each 
parameter direction. The square roots of the diagonal 
elements of the inverse of the Hessian matrix give the 
standard error on each band-power. This procedure re- 
quires only about 30 CPU-mins for ^ 4000 visibilities. 

4.2. Gibbs Sampling Method 

As discussed in iKarakci et al.l ([2013D . the method of 
Gibbs sampling has been applied to interferometric ob- 
servations of the polarized CMB signal in order to recover 
both the input signal and the power spectra. 

The CMB signal is described as a 3np dimensional vec- 
tor, s, of the Fourier transform of the pixelated signal 
maps of rip pixels; s = (..., Tj, Ei, Bi, ...); i = 0, rip — 1. 

The Gibbs sampling method is employed to sample the 
signal, s, and the signal covariance, S = (s s^), from the 
joint distribution P(S, s, dy) by successively sampling 
from the conditional distributions in an iterative fashion 
(Larson et aLlBOOTl: Ikarakci et al.ll2013f) : 



ga+l ^ p(g I (22b) 

After a "burn-in" phase, the stationary distribution of 
the Markov chain is reached and the samples approxi- 
mate to being samples from the joint distribution. 

To determine that the stationary distribution of the 
Markov chain has been reac hed, the Gelman - Rubin 
(GR) statistic is employed (iGelman fc RubinI 119921 : 
iSutter et al]|2012l: IKarakci et al.ll2013t ). For multiple in- 
stances of chains, when the ratio of the variance within 
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each chain to the variance among chains drops to a value 
below a given tolerance, the convergence is said to be at- 
tained. The convergence of the Gibbs sampling is reached 
roughly in 30 CPU-hours. 

5. RESULTS 

5.1. Power Spectra 

The mean posterior power spectra, together with the 
associated uncertainties at each £-bin, obtained by the 
methods of Gibbs Sampling (GS) and Maximum Likeli- 
hood (ML) for the ideal linear experiment, are shown in 
Figure 2. The input power spectra, which are used to 
construct the signal realization, and the spectra of the 
signal realization are also shown in Figure 2. Almost all 
of our estimates fall within 2a of the expected value. 

5.2. Effect of Errors 

In order to estimate a we ran 30 realizations of each 
systematic error simulation for both linear and circular 
experiments. To keep the value of less than 10% 
tolerance limit at r = 0.01, we set the rms values of 
the parameters for gain errors to \grms\ = 0.1, for cou- 
pling errors to \erms\ = 5 x 10""*, for pointing errors to 
6rms ~ 0.7°, for beam shape errors to (rms ~ 0.7°, and 
for cross-polarization errors to firms = 5 x 10~*. 

Figure 3 shows the mean values of a'^^ for beam er- 
rors, averaged over 30 realizations. The results from ML 
and GS methods are in good agreement for both linear 
and circular experiments. In all three cases a^^ ^ 0.1 at 
low i, as expected. Although the cross polarization has 
a much smaller error parameter, its effect on the power 
spectra is comparable to the pointing and shape errors. 
The reason for this is the leakage from TT power into 
BB power that is caused by the off-diagonal elements of 
the beam pattern, whereas the source of a^^ for point- 
ing and shape errors is the EE — BB leakage (Bunn 
2007). 

The mean values of a'^^ for instrumental errors are 
shown in Figure 4. For gain and coupling errors, a"^^ 
is roughly at the 10% level. The main contribution for 
the a comes from the leakage from EE power into BB 
power for gain errors. As in the case of cross-polarization 
errors, despite having a much smaller parameter than 
gain, a^^ ^ 0.1 at low £ for antenna coupling errors 
because of TT BB leakage. 

We simulated the systematics by turning on one error 
at a time. However, in a realistic experiment, all system- 
atic errors act together simultaneously, causing a larger 
effect on the spectra. In order to see this combined effect 
we ran 30 realizations with all the systematic errors dis- 
cussed in previous sections turned on at once. The results 
are also shown in Figure 4. As expected, the combined 
effect is almost twice as large as the individual cases. 

5.3. Comparison to Analytical Estimations 

Analytical estimations for are obtained from the 
quadrature difference of Eq. I15[ normalized by the num- 
ber of baselines. In general a'^^ has a polynomial de- 
pendence on s^. For our interferometer configuration 
is roughly ^ 262. 7/£^. The explicit forms of the un- 
normalized estimations are given in the Appendix. 

In general, our simulated results are larger than the 
estimations in all ^-bins. This is expected because our 



analytical estimations are only first order approximations 
where it is also assumed that the errors associated with 
baselines are uncorrelated, making them lower bounds 
for the estimations. In reality, there is a correlation be- 
tween errors associated with baselines having common 
antennas, a fact that is captured by our simulations. 
Upper bounds for the estimations can be found by un- 
realistically assuming full correlation of errors between 
baselines, where each baseline has the same error. For 
our interferometer design, this corresponds to roughly 65 
times larger values. We expect our results to fall between 
uncorrelated and fully correlated estimations. In order to 
compare our results with the analytical ones, we consider 
the rms values of averaging over the i-h'ms. Figure 
5 shows the ratios of a^^s obtained by ML and GS meth- 
ods to the estimated a^j^^. In most cases, both methods 
are in agreement with the analytical results within a fac- 
tor of 6. 

5.4. Biases in Tensor-to- Scalar Ratio 

The major goal of QUBIC-like experiments is to detect 
the signals of the primordial B-modes, the magnitude 
of which is characterized by the tensor-to-scalar ratio 
r. In this context, it is necessary to propagate the ef- 
fects of systematic errors through to r to assess properly 
the systematic-induced biases in the primordial B-mode 
measurements . 

The shape of the primordial BB power spectrum 
^ I prim insensitive to r but the amplitude is directly 
proportional to r. We can straightforwardly convert the 
amplitude of the systematic-induced false BB into the 
bias in r by writing Cj^p^^^ = r CfJL^ in Eqs. HandUO] 

where CfJ^^^ is the CAMB ([Lewis et al.ll2000[ ) calculated 
primordial BB power spectrum at r = 1. The tensor-to- 
scalar ratios obtained from ideal linear experiment by 
Gibbs sampling and maximum likelihood methods are 
found as res = 0.026 ±0.012 and ruL = 0.006 ± 0.0095, 
respectively. A more conservative estimation for r can 
be obtained without subtracting the lensed spectrum in 
Eq. [9] and by taking only the first bin where the ef- 
fect of lensing is the least; r^g'*^'* = 0.038 ± 0.014 and 

flensed ^ 0.0196 ± 0.011. 

We vary each systematic error individually and also 
consider the cross contributions between each error. In 
realistic observations, all different systematic errors are 
likely to occur at the same time and we need to un- 
derstand their combined effects well. We thus evaluate 
such effects by simulating the systematic errors occurring 
simultaneously during the observation. The individual 
and combined systematic-induced biases in r are illus- 
trated in Figure 6, evaluated by both the GS and ML 
methods based on the simulations performed in the lin- 
ear and circular bases. Both methods demonstrate good 
agreement, within a factor of 2.5. Although the mock 
visibility data are simulated based on only one realiza- 
tion of CMB anisotropy fields, drawn from the power 
spectra with input BB for r — 0.01, the resulting false 
BB band-powers for the different systematic errors are 
expected to be a good approximation for other r values 
since the leading-order false _B-modes are contaminated 
only by the leakage of TT, TE and EE power spectra, 
which are independent of r. 

The simulations show that, due to the leakage of TT 
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Figure 2. Mean posterior power spectra obtained by Gibbs Sampling (GS) for each ^-bin are shown in black. The power spectra 
estimations obtained by Maximum Likelihood (ML) method are shown in blue. Dark and light grey indicate Icr and 2(t uncertainties for 
Gibbs sampling results, respectively. The binned power spectra of the signal realization are shown in pink. Red lines are the input CMB 
power spectra obtained by CAMB for a tensor-to-scalar ratio of r = 0.01. 



signals into BB, even though the cross-polarization and 
coupling errors are very small, e.g., firms = 5 x 10"'' and 
krmsl = 5 X 10~^, the resulting biases in r are compa- 
rable to those induced by relatively larger pointing, gain 
and shape errors. In addition, when increasing the cross- 
polarization and coupling errors by a factor of 10, the 
simulations show that the resulting biases would roughly 
increase by the same factor. As expected, the systematic 
errors are approximately linearly proportional to their 
error parameters. We also find that the combined sys- 
tematic effects (referred to as "c" in Figure 6) would 
increase the biases and their values are consistent with 
the quadrature sum of the individual errors within 10%. 

If we set up an allowable tolerance level of 10% on r, 
where r is assumed to be r = 0.01, for QUBIC-like ex- 
periments the error parameters adopted as in Figure 6 



satisfy this threshold when each systematic error occurs 
alone during observations. But if all the systematic er- 
rors are present at the same time, on average, we require 
roughly 2 times better systematic control on each error 
parameter. Although the tolerance level for r is chosen 
to be ar = 0.1, our results can directly apply to any 
other desired threshold level as long as the linear depen- 
dence of systematic effects on error parameters is a good 
approximation for sufficiently small error parameters. 

6. DISCUSSIONS 

In this work a complete pipeline of simulations is devel- 
oped to diagnose the effects of systematic errors on the 
CMB polarization power spectra obtained by an interfer- 
ometric observation. A realistic, QUBIC-like interferom- 
eter design with systematics that incorporate the effects 
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Figure 3. Beam errors. The values of a-'^^ , averaged over 30 simulations, obtained by both maximum likelihood (ML) method (triangles) 
and the method of Gibbs sampling (GS) (solid dots) are shown. The three rows indicate, from top to bottom, pointing errors with 
Srms ~ 0.7°, beam shape errors with Crms ~ 0.7°, and beam cross-polarization with ^irma = 5 X 10~*. Left panel shows (red) 
and (blue). Middle panel shows (red) and (blue). Right panel shows (red) and (blue). Linear and circular 

experiments are shown by solid and dashed lines, respectively. 



of sky-rotation is simulated. The mock data sets are an- 
alyzed by both the maximum likelihood method and the 
method of Gibbs sampling. The results from both meth- 
ods are found to be consistent with each other, as well 
as with the analytical estimations within a factor of 6. 

In order to assess the level at which systematic effects 
must be controlled, a tolerance level of = 0.1 is chosen. 
This ensures that the instrument is sensitive enough to 
detect the B-signal at r = 0.01 level (O'Dea et al. 2007). 
We see that, for a QUBIC-like experiment, the contami- 
nation of the tensor-to-scalar ratio at r = 0.01 does not 
exceed the 10% tolerance level in the multipole range 
28 < i < 384 when the Gaussian-distributed system- 
atic errors are controlled with precisions of |(7rms| = 0.1 
for antenna gain, |ei.ms| = 5 x 10^** for antenna cou- 
pling, Srms ~ 0.7° for pointing, (rms ~ 0.7° for beam 
shape, and ^rms = 5 x 10""* for beam cross-polarization 
when each error acts individually. However, in a real- 
istic experiment all the systematic errors are simultane- 



ously present, in which case the tolerance parameter of 
r roughly reaches the 20% level, suggesting that better 
control of systematics would be needed. 

Apart from the systematics presented in the paper, we 
also ran simulations to analyze the effects of uncertainties 
in the positions of the antennas. In order to have an 
effect on the order of a^^ = 0.1, we found that the 
uncertainty in the position of each antenna should be on 
the order of 50% of the length of the ww-plane. Since 
such an error is unrealistically large, we conclude that 
the effect of antenna position errors on power spectra is 
negligible in an interferometric observation. 

We have shown that a QUBIC-like experiment has 
fairly manageable systematics, which is essential for the 
detection of primordial i?-modes. Since our interferome- 
ter design has a large number of redundant baselines (ap- 
proximately 10 baselines per visibility), as a further im- 
provement, a self-calibration technique can be employed 
to significantly reduce the level of instrumental errors 
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Figure 4. Instrumental and combined systematic errors. The values of a , averaged over 30 simulations, obtained by both maximum 
likelihood (ML) method (triangles) and the method of Gibbs sampling (GS) (solid dots) are shown. Top row: antenna gain with \grmB | = 0.1. 
Middle row: antenna couplings with |erms| = 5 X lO"**. Bottom row: combined effect of beam and instrumental systematic errors. Left 
panel shows (red) and (blue). Middle panel shows (red) and (blue). Right panel shows a-^^ (red) and (blue). 

Linear and circular experiments are shown by solid and dashed lines, respectively. 



(ILiu et al.ll2010l) . 
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APPENDIX 



Following IBunnl (|2007D . we obtain first order approximations for the AC;?^ given, for a single baseline, in Eq. 1151 
For a baseline lying on the x-axis, the matrices in Eq. [T2] are given as 



TT 




N 



EE 



TE 




N 



TB 




The covariance matrix can be written in block- matrix form as A^t 



Mo Ml 
rt 



Ml M2 

c 



where 



Mo 



c 



TB-, 



(JEE^2 



(JBB^2 



For a baseline lying in an arbitrary direction, these matrices must be transformed as Mq R ^MqR and Nxy — ) 
R~"'^NxyR, where R is the rotation matrix given in Eq. 1141 The resulting expression will, then, be averaged over Q. 



For instrumental errors (5v = E • v, which gives Mi Mq • and M2 = E • Mo • E^. 



GAIN ERRORS 



Linear Basis 



91 ^1(91+9-2+ 91* + 9i*), 92 -1(91-92+ " 911 



71 = 1(9? + 9Y), 72 = 1(9? - 9T), 73 = ^(ff^^ + 9^), 

/ 71 73 \ 

Ef" = 71 + 72 

linear \ I ' I I 

V 71 - 72 / 
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(ACJ^J^ ^ (6i?e{7i}2 + |i?e{72}^ - i/m{72}^)(C^^)2 + (2^71^ + ^172^)^^^^^^ 

{AC,^Zr = (8i?e{7i}2 + 4i?e{72}^)(C^^)^ 

(Aa^„^j2=F|72|2(t7^T 

(A(7J^J2 ^ 1 (3ije{72}2 - /m{72}^)(C^^)^ + i(l72|' + 8Re{^,y)C^^C^^ 
{AC^Zf = Re{^,Y{C^^f + (|72p + 2Re{^,Y)C^^C^^ 

Circular Basis 

(91 ON 

^lo. = 51 i52 

V -192 91 J 



{l^Cfl,f = me{9^Y{C^^f 

i^c^zr = {iim{92r - hRe{92y){c^^r + 'm'c^^c 



COUPLING ERRORS 



Linear Basis 



ei = ^(el + e| + ef + ), = ^{e\ - - e{* + ) 



ei 

9 _ ( ^. 

"'imear 



EJ:°ZT= I ei+e2 -e3-e4 



\ ei - £2 £3 - €4 / 
(AC',^X)2 = (3i?e{ei}2 - 7m{ei}2)(C^^)2 + |ei|2(7^r^£i? 

(A(7f^j2 ^ 2(|ei|2 + |e2|2)C^^C^^ 

(A(7,^^J2 ^ 2(|ei|2 + |e2|2)C^r^Bi3 + 2^(|ei|2 + |e2nC^^C^^ 

{AC^Z? = ideiP + |e2p)C^^C^^ + i(3i?e{ei}2 + 3i?e{e2}2 - /m{ei}2 - 7m{e2}2)(C^^)2 



Circular Basis 



ei ie2 



j,co«pims = I e 

n.mi.l.n.r I J- 



ie2 



POINTING ERRORS 

Defining Sr^ as the deviation of the k^^ antenna's pointing center, we can write, to the first order, 

A^{r)A^{r) = exp[-{r - (7Sf/2a% 

where S = {Srj + Srk)/2a. 

SVz = -ia j d2kZ(k)[i2(k - 27ru)]*[(k - 27ru) • 8z\ 

{Vx5V^) = and {5Vx5V^) = • 5y) {VxV*) . 
Linear Basis 

Si = -j^{SQ + 6u), S2 = -{Sq-Su) 

(ACJ^J^ ^ i|5^|2(cTi^)2 + i(4|5^|2 + \s,nc^^C^^ 
(Aa^™^j2 = (l'5iP + ^l^2p)(C^^)^ 

Circular Basis ^2 = 
SHAPE ERRORS 

The product of two elUptic Gaussian beams can be written as a single eUiptic Gaussian: 

A^(v)A''(r) = exv [_ (^£^£A±1£^^ _ (ycos/3 - a;sin/3)^ ' 
^ ^ ^ ' ^[ 2{a + (T^)2 2(<T + ay)^ J ' 

where /3 is the angle between the major axis of the resulting ellipse and the x-axis. 

5Vz = / cf^Z{k)[{^Az){k - 27ru)]* 

where Az{x,y) = a;2(Cf cos^ /3 + sin^ /3) + ^^(g^ cos^ /3 + Cf sin^ /3) + a;y(Cf - C^)sm2p, and 
The only non-vanishing integrals in the covariance matrix are: 
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J |i2|2 = J A^x^A-^y = J A^y^A^)* = ^Tra\ J {x^A^l"" = J \y-^A^\^ = ^7^c7^ 

j (^)(^)* = j |iyX2|2 ^ 1^^6_ 

Averaging over /3 we get (Vx^V^^) = -Cf {VxV*) and (^Vx^l^^) = (2Ci^Cr + C2CI) (VxV*) . 
Linear Basis 

(A4^JJ2 = (10C?+ + 2C|+)(C^^)2 

(ACji^j2 ^ (7^2^ + 3^2_ + + (3C?+ + K?- + CI+ + iCl_)C^^C^^ 

(Aa^„^j2 ^ (ioC?+ + 5CL + 2CI+ + C|_)(C^^)2 

(Aa^^j2 ^ (2^2_ + c|_)^(C^^)2 + (12C?+ + 12C?_ + 4CI+ + a_)-^C^^C^^ + (10C?+ + 5C?_ + 2C|+ + Cl_)(C^^)2 

(AC-J^J^ = 3^2_(cTi5)2 + 1 (2C2_ + ^2_^C^TCEE + (3^2^ + 1^2_ + ^2^ + 1 ^2_ )^C'^T 
(ACf„Sj2 = (5^2_ + 1^2_ +^(3^2^ +^2^))(ci^i^)2 

Circular Basis Ci- = 
CROSS-POLARIZATION 

The only non-vanishing integrals in the covariance matrix are: 

J |i2|2 = na^, J |A2^^2^2(/.|2 = J |A2^^2(/.|2 = 7ra^ 
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\ ^ rms ) 


= 2/i?C^^C^^ 












= 2(M? + Mi)C^^C'^^ 
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Circular Basis 





1 9 9 9 
P P P P 

61 = u+^Q sm2(p, 6Q — Li+^I sm2d) + iu-^U cos24i, 6U = —iii-^Q cos2(j) 
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( ,,2 nTT o.,2 r~iEE\r<EE 
= — ZIJ,_U jO 


(AC^g)^ 













